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Abstract 



In the past, several efficient methods have been developed to solve the 
Schrodinger equation for four-nucleon bound states accurately. These are 
the Faddeev-Yakubovsky, the coupled-rearrangement-channel Gaussian-basis 
variational, the stochastic variational, the hyper spherical variational, the 
Green's function Monte Carlo, the no-core shell model and the effective in- 
teraction hyperspherical harmonic methods. In this article we compare the 
energy eigenvalue results and some wave function properties using the realis- 
tic AV8' NN interaction. The results of all schemes agree very well showing 
the high accuracy of our present ability to calculate the four-nucleon bound 
state. 

21.45.-Fv, 24.10.-i, 27.10.-Fh 
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Typeset using REVT^ 



I. INTRODUCTION 



Recent advances in computational facilities, together with the development of new meth- 
ods and refinements upon older ones, allow very precise calculations for few-body systems. 
These advances are especially remarkable in nuclear physics considering the complexity of 
the nuclear interaction. The three-nucleon(3N) bound-state and scattering-state ||^-|^ 
problems are rigorously solved using realistic nuclear potentials These calculational 

schemes are mostly based on a partial wave decomposition. Stochastic and Monte Carlo 
methods for bound states, however, are performed directly using position vectors in configu- 
ration space. Also in momentum space the first steps have been taken to avoid partial wave 
decompositions in both two-nucleon (NN) |]10| and 3N |lTT| , p^ systems. Benchmark calcula- 



tions based on different algorithms for the 3N continuum both below [0,|T1] and above 
the deuteron threshold already exist. 

The complicated calculation of few-body continuum states can be avoided in the eval- 
uation of reaction cross sections, even in the presence of realistic forces In fact the 
transition strength can be calculated in an alternative way, where only bound state tech- 
niques are needed . 

There are a few analytical solutions of 3N bound states |jl8| for square-well potentials. 



against which numerical solutions have been checked, but they are far from possessing the 
complexity of realistic nuclear forces. In the four-body system we are only aware of bench- 
mark calculations for four bosons [Q. 

Benchmark calculations are extremely useful to test methods as well as calculational 
schemes. They are also often of interest for a general readership, since they may help to 
solve analogous problems in other fields. We think that this is particularly the case for the 
quite complex four-fermion system. Here we would like to address the four-nucleon (4N) 
bound-state problem using the AV8' NN potential which is a simplified, reprojected 
version of the fully realistic Argonne AV18 model P], but still has most of its complexity, 
e.g., the tensor force is built in. 

In Sec. 2 the different methods are briefly introduced and the results are presented in 
Sec. 3 together with a brief summary. 

II. METHODS 

In order to solve the bound four-nucleon system we employ the Faddeev-Yakubovsky 
equations(FY) [p0|-|23[], the Coupled- Rearrangement- Channel Gaussian-Basis Variational 
Method (CRCGV) |2§-|ll, the Stochastic Variational Method(SVM) with correlated Gaus- 
sians p2|-p5[1, the Hyperspherical Harmonic Variational Method(HH) p6|-p[1, the Green's 
function Monte Carlo(GFMC) ||,g3|,|l9|,|4l method, the No-Core Shell Model(NCSM) 



1^-^, and the effective interaction hyperspherical harmonic method(EIHH) [Q. The 



various procedures are briefly described below. 



A. Faddeev-Yakubovsky Equations 



The 4N Faddeev-Yakubovsky equations in momentum space are [21-23 
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^2=Gotl2P[(l-P34)V^l+^2], 



(1) 

(2) 



where ipi and 1/J2 are Yakubovsky components and ti2 is the two nucleon transition matrix 
determined by a two nucleon Lippmann - Schwinger equation. P, P34 and P are permutation 
operators: P = P12-P23 + -P13-P23, P = -P13-P24, where the Pij are transpositions of particles i 
and j. The fully anti-symmetrized wave function is 

^ = [1 - (1 + P)Pu]{l + P)^Pi + (1 + P){1 + P)i^2. (3) 

The Yakubovsky equations are decomposed into partial waves. We truncate the partial 
waves at a 2-body total angular momentum j < 6, all other orbital angular momenta at 
Zj < 8 , and the sum of all angular momenta at X)j h < 12. In this truncation we keep 1572 
angular momentum and isospin combinations (often called channels). This is sufficient to 
guarantee convergence of our results as given in Sec. 3. The diagonalization is performed by 



a modified Lanczos method |49]. Recent results for realistic NN potentials, and including 



three-nucleon forces, are given in [E^]. 



B. Coupled-Rearrangement-Channel Gaussian-Basis Variational Method 

The coupled-rearrangement-channel Gaussian basis variational method was proposed by 
Kamimura |^ to solve the Coulombic three-body problem of the muonic molecular ion 



((it/i)"*", within an accuracy of seven significant figures for the energy of the very loosely 
bound J = V = 1 state; this accuracy was required for the comparison with experimental 
data on the muon catalyzed fusion cycle. Use of basis functions that spanned all the three 
rearrangement Jacobian coordinates was essential to the high-precision calculation. The 



method was also applied to three-nucleon bound states p5| , p6| and was found to accomplish 
a much more rapid convergence in the binding energy with respect to the number of the 
three-body angular momentum channels (see Fig. 5 of p^). 

The method was also successfully used to make another high- precision Coulomb three- 
body calculation of the anti-protonic helium atom [p + He^"^ + 6~) in highly excited meta- 



stable states with J 35 The calculation agreed with the high- resolution laser spec- 
troscopic data within seven significant figures so that the mass of antiproton was derived 
to two orders of magnitude better precision than published values. The method has been 
useful in four-body calculations of the structure of light hypernuclei with realistic YN and 
NN interactions 

The total four-body wave function is described as the sum of amplitudes of the 
rearrangement- Jacobian-coordinate channels with the LS coupling scheme : 

^.M = E^f^e^+E^f^*i^^ (4) 

cx a 

where anti-symmetrized basis functions are described with quantum numbers a = 
{nl, NL, A, uX, I, ss'S, t} by 



^L^^ = ^{[[[0„Kr)fe(R)]A Vux{p)]j 
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(5) 



^i''^=A{[[[Mr')^NLm]AV.x{p% 



(6) 



We employ K-tjpe coordinates r = xi — X2, R = (xi +X2)/2 — X3, p = (xi +X2 + X3)/3 — X4 
and H-tjpe ones r' = xi — X2, R' = X3 — X4, p' = (xi + X2)/2 — (X3 + X4)/2. A is the four- 
nucleon antisymmetrizer and x's and ?7's are the spin and isospin functions, respectively. 
The functional form of 0„;(r) is taken as 



Pnlm 



;r)=r'e-('-/-)'l^„(f) 



(7) 



where the Gaussian range parameters are chosen to lie in a geometrical progression 
(r„ = ria"~^;n = 1 ~ 'n.max); and similarly for the other functions ip and if. This man- 
ner of choosing the range parameters is very suitable for describing both the short-range 
correlations and the long-range asymptotic behavior precisely p5|j31[] . 

Eigen-energies and wave-function coefficients C's are determined by solving the 
Schrodinger equation with the Rayleigh-Ritz variational principle. It is to be emphasized 
that truncation is not made for the partial waves of the NN interaction, in contrast to the 
Faddeev-Yakubovsky method, but is done only for the angular momenta of basis functions, 
as in most variational methods. This makes it possible to accomplish a very quick conver- 
gence; the result in Sec. 3. uses l,L,X < 2 (this is the same as in the case of the three 
nucleon bound states, mentioned above). For instance, this amounts to 100 channels for the 
calculation in Sec. 3. 



C. Stochastic Variational Method 



The Correlated Gaussian trial function is written in the following form [p2l-|3^: 



(8) 



where A is the antisymmetrizer, x stands for a set of A — 1 intrinsic coordinates 
(xi,X2, . . . ,xa--i) and ^sMi {^taIti) is the spin (isospin) function of the A-particle system. 
These functions are constructed by successively coupling the spin (isospin) of the nucleons 



e5Mi=[[[Xi(l)Xi(2)].,,Xi(3)] 



S123 • • -Is Mi 



(9) 



(similarly for the isospin part). The non-spherical (orbital) part of the trial function is 
represented by a successively coupled product of spherical harmonics 
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LMiW 



[F,,(X1)1^,(X2)], n,(x3) 



LMi 



(10) 
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The index i in the above equation stands for a label to distinguish the different possible 
intermediate coupling schemes as well as the different possible total spin and total orbital 
angular momentum. 

The expansion over the partial waves has to be truncated, but the correlations included 
in the Gaussian part, exp(— ^xAjx), make the trial function flexible enough, so these trunca- 
tions are not expected to seriously affect the accuracy. In the present calculation we included 
all partial waves up to 

j:h<4. (11) 

i=l 

The trial function contains A{A--l)/2 nonlinear variational parameters. The total spin, total 
orbital angular momentum and intermediate coupling quantum numbers are also variational 
parameters in the sense that one has to include all possibilities which improve the energy. 
We have a large number of parameters to be optimized and it is not at all clear how to select 
the optimal quantum numbers. 

This variational basis is nonorthogonal, none of the components is indispensable, and 
one can replace a component by a linear combination of others. This gives us an excellent 
opportunity to use a stochastic optimization procedure. To optimize the variational basis 
we used the "stochastic variational method" |]32|-p5|. In the SVM one searches for the best 
wave function by a random trial and error procedure. Random trial functions are generated 
and their energies are compared. Randomness in this case means that the quantum numbers 
and the nonlinear parameters are random numbers. Trial functions giving the lowest energy 
are selected as basis states. Details and various applications of the approach can be found 



m 



The number of basis states used in the calculations is about 150 for the triton and 300 
for the alpha particle. Very small bases already give quite acceptable results, for the alpha 
particle, for example, 50 basis states give the binding energy within 1 MeV. The SVM results 
seem to be convergent in the model space defined with 300 basis states and the partial-wave 
truncation with J2t=i k < 4. We have tried to increase the accuracy by adding 700 more 
states and by including the partial waves up to J2i=i^ k ^ Q but the results are practically 
unchanged. The 1000 bases give only 2 keV gain in energy. The enlargement of the basis 
improves the expectation values, especially that of the kinetic energy operator, but this 
change is canceled by a similar change in the central potential. We think that the upper 
bound provided by the SVM calculation is very close to the exact energy. The accuracy 
achieved with few basis dimension is due to the use of the correlated Gaussian basis and the 
efficient optimization procedure. 



D. Hyperspherical Harmonic Variational Method 

The Hyperspherical Harmonic (HH) functions constitute a general basis for expanding 
the wave functions of an A-hodj system [p6|-p^ . Very precise results can be obtained for the 
three-nucleon bound state ||39|. In the HH variational method, the wave function is written 

as 

^ = j:mpH''^. (12) 
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where p is the hyper-radius. The quantities are fully anti-symmetrized HH-spin-isospin 
functions of quantum numbers p = {n,m, ii, £2, ^3, ^12, L, s, s' , S,t,t' ,T} constructed using 
the ii'-type Jacobi coordinates xi^2,3- Explicitly, 



(sin /5)2-p:'''+^ (pQg 2/3)P™^+^'''''^ (cos 27) 



Vtil2HS) MA) 



X 



[Ye,{x,)Ye,{x2)],^Ye,{xs)] L(12)x(3)l y(4) 



s) J J:, 



(13) 



where cos/3 = x^/p, cos 7 = X2/ {psin P) and x and r] denote spin and isospin functions, 
respectively. Moreover, z/ = £1 + £2 + 2m + 2 and P^'^ are Jacobi polynomials (the integers n 
and m range from zero to infinity). The coefficients M^(p) depend on the hyper-radius p = 

\J {xiY + (0:2)^ + (3:3)^ and can be determined by solving a set of second-order differential 
equations derived from the Rayleigh-Ritz variational principle. For A = 4 the necessary 
matrix elements of the potential have been calculated by exploiting the techniques discussed 
in Ref. ||. 



The main difficulty in applying the HH technique to nuclear systems is the very slow 
convergence of the expansion due to the strong repulsion between the particles at short 
distances. In the A = 4 case, it has been found convenient to separate the HH states 
in different classes and to study the convergence by including the states of one class at 
time The adopted criterion has been to first include the HH functions describing two- 
body correlations and, successively, those incorporating three- and four-body correlations. 
Moreover, the HH functions having the lowest orbital angular momentum quantum numbers 
ii {i = 1, 2, 3) have been included first. 



E. Green's Function Monte Carlo Method 

Green's function Monte Carlo methods use stochastic sampling to evaluate path integrals 
of the form: 

^0 = lim ^(r) , (14) 
^(r) = e-(^-^o)"^T , (15) 

, (16) 



where is an approximate trial function obtained in a variational or in an approximate 
constrained-path GFMC calculation and we have introduced a small time step, r = nAr. 
An approximate expression for the propagator, 

G(R, R') = (R|e-(^-^«)^nR') , (17) 

with error proportionate to at least the second power of At is used. In the present work we 
use a symmetrized product of exact two-body propagators, which has error proportionate 
to (Ar)3 igg. 

Green's function Monte Carlo calculations for light nuclei with spin-isospin dependent 
interactions sample the particle coordinates while explicitly summing over the spin-isospin 
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degrees of freedom |H2| . The first alpha particle calculations including L ■ S terms employed 



the Reid V8 interaction ^3[. The chief advantage of these methods is that they can be 
extended to larger nuclei. More computationally efficient versions of the algorithm have 
been introduced and calculations extended up to A=8 [^,0. 

Convergence of the ground-state energy is governed by the spectra of the Hamiltonian. 
Calculations reported here were performed to r = 0.12 MeV~^. Since the first excited state 
of ^He is above 20 MeV, any errors in are damped out by at least exp(— 2.4), an order 
of magnitude. In fact our studies show that the errors in correspond to much higher 
excitation energies and {H) converges by r ~ 0.05MeV~^. 

The GFMC method allows us to compute mixed expectation values of the form 
{^{t)\0\^t)- For H, this gives the exact ground state energy if r is large enough. Ex- 
pectation values of other quantities, such as pieces of the Hamiltonian, are often obtained 
through a linear extrapolation in the error of the trial wave function: 

(^(r)|C|^(r)) ^ 2(^(r)|C|^T) - (^t|C|^t), (18) 
though it is possible to go beyond this approximation. 



F. No-Core Shell Model Method 



The NCSM is an approach applicable to both few-nucleon systems as well as to light 
nuclei |^5| . The calculations are performed in a finite model space in the harmonic-oscillator 
(HO) basis. The model space (P) is spanned by states with the total number of HO quanta 
N < A^max- The Hamiltonian, 



H=T+V , 

is modified by a HO center-of-mass potential. Thus, we work with 



A 

= E 

1=1 

A 

+ E 

i<j=i 



Pi I l ri2 2 

— mi 2 r,- 

2m 2 



(19) 

(20) 
(21) 



As the NN potential depends on the relative coordinates, the added HO term has no infiuence 
on the internal motion in the full space. The effective Hamiltonian, appropriate to the finite 
P-space, is derived by the hermitian version of the Lee-Suzuki method |4^. In general, the 
effective Hamiltonian is an A-body operator. We make an approximation by using just a 
two(three)-body effective interaction, which is obtained by applying the Lee-Suzuki approach 
to the two(three)-nucleon system using with the sums restricted to two(three) nucleons, 
but with the A in the interaction term kept fixed to, e.g., A = 4 for ^He. Consequently, we 
deal with a two(three)-nucleon system bound in a HO potential. The effective two(three)- 
body interaction, then replaces the interaction term in H^. We note that the effective 
interaction, by construction, converges to the original bare interaction as the basis space 
is increased and, thus, the NCSM calculation converges to the exact solution with the 
basis-space enlargement. In fact, it converges much faster than the corresponding bare 



8 



interaction calculation performed in the same basis. Eventually, the A-nucleon P-space 
calculation can be performed either in a Slater-determinant single-particle HO basis or in 
a properly antisymmetrized Jacobi-coordinate HO basis. The latter is used in the present 
^He calculation. In the past, we applied this approach successfully to the ^He interacting by 
the CD-Bonn NN potential. It turns out that the convergence with the AV8' is significantly 
slower. The limitation to a two-body effective interaction is inadequate in the P-spaces that 
we could access (A^max = 18). Therefore, we performed the calculations using the three-body 
effective interaction. The mean values of the different operators were calculated using the 
corresponding effective operators computed within the Lee-Suzuki approach in a two-body 



approximation using the formula derived in Ref. E7 



G. Effective Interaction Hyperspherical Harmonic Method 

Similarly to the preceeding method the EIHH approach introduces a two-body effective 
interaction Veff The division of the total HH space in P and Q spaces is realized 

via the HH quantum number K (P(Q) space: K < {>)Kmax)- Two powerful algorithms 
recently developed for the construction of symmetrized HH functions are employed |]50| , |5T 
In hyperspherical coordinates the total Hamiltonian is written as 



where p is the hyper-radius and Ap contains derivatives with respect to p only. The grand- 
angular momentum operator is a function of the variables of particles A and (A — l) and 
of Ka-2 the grand angular momentum operator of the {A — 2) residual system ^2|. Then 



from the total Hamiltonian one can extract a "two-body" Hamiltonian of particles A and 
{A -I) 

^2(p) = ^^ + V^A(A-i), (23) 

which, however, contains the hyperspherical part of the total kinetic energy. Since the HH 
functions of the {A — 2) system are eigenfunctions of K\_2 one has an explicit dependence 
of H2 on the quantum number Ka-2 of the residual system, i.e. H2 H2^'^. Applying 



the hermitian version of the Lee-Suzuki method |46] to H2 one gets an effective Hamiltonian 



H2eff- The effective interaction V^jf is obtained from 



K/r(p)=<r(p)-i^ (24) 



2m p^ 

This Veff replaces Vij in Eq. (|2^ when we project the solution on the P-space. This effective 
potential has the following property: Veff — > Vij for P — > 1. Due to the "effectiveness" of the 
operator the solution of the Schrodinger equation converges faster to the true one. The HH 
formulation leads to various advantages: (i) Veff itself is p dependent, therefore it contains 
some information on the "medium" , (ii) because of the above mentioned Ka_2 dependence 
the (A-2) residual system is not a pure spectator, and (iii) an additional confining potential 
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is not needed, since the presence of p in Eq. ( p3D automatically confines the two-body 
system to the range < r^_(A-i) < "\/2p. We would like to point out that K//(-^max) can 
be viewed as a kind of momentum expansion, since the short range resolution is increased 
with growing Kmax- As discussed for the NCSM approach one obtains a better convergence 
for the calculation of mean values introducing corresponding effective operators. Of course 
for the calculation of the mean value of the Hamiltonian, i. e. Ei,, one already makes use of 
an effective operator, namely H^/ff^- 



III. RESULTS 

The AV8' interaction appears to be an ideal test potential to compare the different 
calculational schemes. It is derived from the realistic AV18 interaction |^ by neglecting the 
charge dependence and the terms proportional to and {L ■ S)"^. Furthermore, in this work 
we omit the electromagnetic part of the interaction. The potential is local and its spin and 
isospin dependences are represented by operators. Because of its form it is tractable for all 
of the calculational schemes described above. 

The potential consists of 8 parts: 

V{r) = K(r) + K(r)(r ■ r) + V^{r){a ■ a) + Kr(r)(a ■ a){r ■ r) 
+Vt{r)Si2 + Vtrir)Su (r • r) 

+n(r)(L ■ S) + Vbr{r){L ■ S) {r ■ r) (25) 
= j:V,{r)0,, (26) 

i=l 

where (a ■ a), (r ■ t),Si2 and (L ■ 5*) stand for spin-spin, isospin-isospin, tensor and spin- 
orbit interactions 0] respectively, and Vi{r) are radial functions of Yukawa- and Wood-Saxon 
types. The AV18 and AV8' are defined with ^^/m7v=41.47108MeV /m^, computed from the 
average of the proton and neutron masses. Most of the results reported here were obtained 
using the traditional value of 41.47; this results in a change in (H) of only ^ 2.6 keV, far 
less than the estimated errors in the various methods. 

First, we compare the binding energy results Eh, the expectation values of the kinetic 
and potential energy and the radii in Table |. We find good agreement for E^ within 3 digits 
or within 0.5 %. This is quite remarkable in view of the very different techniques and the 
complexity of the nuclear force chosen. Except for NCSM and EIHH, the expectation values 
of T and V also agree within 3 digits. The NCSM results are, however, still within 1 % and 
EIHH within 1.5 %of the others, but note that the EIHH results for T and V are obtained 
with bare operators. The uncertainty in the NCSM results is of the same size, i.e., IMeV, 
as that for the GFMC. Finally, the given radii are also in very good agreement. 

The HH calculation includes about 4500 states with £ = £i -|- £2 + ^3 < 6. The states 
with C = 6 give a contribution to the binding energy of approximately 0.04 MeV. It is to 
be noticed that the HH spin-isospin states ^^^^ having C < 6 but constructed with the 
H-type Jacobi coordinates are linearly dependent on those considered in the expansion and 
therefore it is unnecessary to include them. The contribution of ^^^^ (and ^^^^) to the 
binding energy with £ > 8 has been estimated to be approximately 0.01 MeV. 
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The errors quoted for the GFMC results are just the Monte Carlo statistical errors. 
Various tests show that the energy is converged to at least this accuracy for changes in Ar 
or the maximum r. There should be no other sources of systematic error in this simple test 
case. 

The NCSM binding energy result is based on extrapolation from calculations using the 
three-body effective interaction in model spaces up to A'^max = 16 in the HO frequency 
range hQ = 16 — 43 MeV. The mean values of different operators, evaluated for A^max = 16 
consisting of 2775 basis states and hfl = 28 MeV, were computed using effective operators 
as the use of bare operators is completely insufficient, in particular for the Vc{r) and T. Note 
that we have here {T^s) + (Ks) close, but not exactly equal to (Hcs), due to approximations 
used. Overall, the NCSM results are less accurate than the other methods. The NCSM 
convergence rate is rather slow for the AV8'. However, the method is flexible to handle also 
non-local realistic potentials like the CD-Bonn with a faster convergence rate due to a softer 
repulsive core. The advantage of the method is its applicability to the p-shell nuclei. 

The EIHH calculation is carried out with K^ax = 20, (about 3000 HH states). The 
error estimate is based on the convergence with respect to K^ax, i-e. difference of results 
for Kmax = 18 and 20. An inspection of Table I shows that Eb and radius are converged to 
a very high precision [Ef,: 0.04 %, radius: 0.007 %, not shown in Table I ). On the contrary 
(T) and {V) still change by about 1 % from K^ax = 18 to K^ax = 20. Of course, by 
construction of the EIHH method, also (T) and {V) have to converge to the true result. 
In order to have a higher precision one can proceed in two ways: (i) increase of Kmax, (ii) 
use of effective operators. Particularly advantageous is the use of effective operators, since 
it allows us to make rather precise calculations with a small number of basis functions (see 
discussion of EIHH result for Fig. 1). As Table I shows it is not necessary to use effective 
operators for long-range observables like the radius, while observables that contain short 
range information (high momentum contributions), like (T) and {V), should, in principle, 
be calculated with effective operators. 

A more detailed test of the wave function is to evaluate the expectation values of the 
eight individual potential energy operators in eq. (^). The results are shown in Table 
|T[ The agreement is, in general, rather good and well within 1 %, except for NCSM with 
discrepancies up to 6 % but they are generally 4% or less. In the case of the CRCGV, the 
expectation values for the spin-orbit operators are a bit off from the rest, but again still 
within 4 %. There are no results given for the EIHH. 

Table |TT| shows the expectation values of the sum of the first 4 operators in eq. ( pSj ) 
(called Central), of the two tensor operators and of the two spin-orbit operators. Again, no 
results are given for the EIHH. Except for the NCSM with differences up to 3.2 %, all the 
values agree quite well each other. 

As a further property of the wave function we consider the NN correlation function 

C(r) = (^|5(f-n2)|^), (27) 

where fi2 = ri — ^2. It is apparently normalized as Att J C{r)r'^dr = 1. The results for the 
various calculational schemes, except for the GFMC are shown in Fig. 1. The agreement 
among the FY, CRCGV, SVM, HH and NCSM is essentially perfect. For the EIHH it is 
necessary to use an effective operator in order to obtain good convergence also for r < 1.2 
fm. Due to the use of rather unsophisticated computers, the EIHH calculation for C{r) is 
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performed with the rather low K^ax value of 12 (about 400 HH states); however, a rather 
good agreement with the other methods is already obtained at this low value. 

Finally, we show in Table ^ the probabilities for finding the three different total orbital 
angular momenta in our 4N model system. The agreement among the different methods is 
very good with a small excursion in NCSM. 

To summarize, we have demonstrated that the Schrodinger equation for a four-nucleon 
system can be handled quite reliably by different methods leading to very good agreement in 
the binding energy, in expectation values of the kinetic and potential energies and in simple 
wave function properties. The AV8' NN potential encompasses most of the complexity of 
realistic NN forces and, thus, the benchmark calculations are highly nontrivial and demon- 
strate the maturity and reliability of various methods. These results are good foundations 
for further investigations of nuclear structure for more complex systems and/or for other 
NN interaction models. 

We have chosen the AV8' potential because it can be handled without any approximation 
by all of our methods. More realistic NN potentials such as AV18 p[, CD-Bonn |^ and 
Nijmegen I, II ||^ pose additional difficulties for at least some of the methods. There are new 
operator forms with higher order derivatives or very strong nonlocalities. Also some of the 
potentials are defined partial wave by partial wave. 

Whereas in the four- body system the FY and NCSM schemes can handle all types 
of NN potentials directly, the GFMC method relies on AV8' and treats the difference to 
AV18 in perturbation theory. The SVM can in principle treat any local potential, such as 
AV18 , but the {L ■ Sf terms require additional computational effort. Also the remaining 
methods, HH , CRCGV and EIHH can handle more complicated potentials, although at 
present applications have been restricted to local potentials. GFMC, NCSM, SVM and 
EIHH have already obtained solutions for A > 4, whereas FY up to now has been restricted 
to y4 < 4. An advantage of the methods, CRCGV and EIHH, is that they do not need as 
heavy computational facilities as the other methods. 
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TABLES 



TABLE L The expectation values (T) and {V) of kinetic and potential energies, the binding 



energies Ej, in MeV and the radius in fm. 



Method 


(T) 


{V) 




W) 


FY 


102.39(5) 


-128.33(10) 


-25.94(5) 


1.485(3) 


CRCGV 


102.30 


-128.20 




-25.90 




1.482 


SVM 


102.35 


-128.27 




-25.92 




1.486 


HH 


102.44 


-128.34 




-25.90(1) 


1.483 


GFMC 


102.3(1.0) 


-128.25(1.0) 




-25.93(2) 


1.490(5) 


NCSM 


103.35 


-129.45 




-25.80(20) 


1.485 


EIHH 


100.8(9) 


-126.7(9) 




-25.944(10) 


1.486 


TABLE IL Expectation values of the eight potential operators in Eq.(|25|) in MeV 


Method 


{Vc) 


(K) 






FY 


16.54 


-5.038 




-9.217 




-57.55 


CRCGV 


16.54 


-5.035 




-9.215 




-57.51 


SVM 


16.54 


-5.036 




-9.213 




-57.51 


HH 


16.57 


-5.034 




-9.255 




-57.59 


GFMC 


16.5(5) 


-5.03(6) 




-9.21(7) 




-57.3(5) 


NCSM 


16.16 


-4.92 




-9.77 




-57.89 


Method 


iVt) 


(Vtr) 


m 


(Vbr) 


FY 


0.707 


-69.06 




10.79 




-15.50 


CRCGV 


0.708 


-68.99 




10.60 




-15.30 


SVM 


0.707 


-69.03 




10.78 




-15.49 


HH 


0.702 


-69.03 




10.76 




-15.46 


GFMC 


0.71(3) 


-68.8(5) 




10.62(15) 




-15.40(15) 


NCSM 


0.68 


-69.13 




11.23 




-15.80 



TABLE HL Expectation values of potential energy operators in MeV. 



Method 


Central 


Tensor 


Spin-orbital 


FY 


-55.26 


-68.35 


-4.72 


CRCGV 


-55.22 


-68.28 


-4.70 


SVM 


-55.23 


-68.32 


-4.71 


HH 


-55.31 


-68.32 


-4.71 


GFMC 


-55.05(70) 


-68.05(70) 


-4.75(5) 


NCSM 


-56.43 


-68.45 


-4.57 
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TABLE IV. Probabilities of total orbital angular momentum components in [%] 



Method 


S-wave 


P-wave 


D-wave 


FY 


85.71 


0.38 


13.91 


CRCGV 


85.73 


0.37 


13.90 


SVM 


85.72 


0.368 


13.91 


HH 


85.72 


0.369 


13.91 


NCSM 


86.73 


0.29 


12.98 


EIHH 


85.73(2) 


0.370(1) 


13.89(1) 
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FIGURES 

FIG. 1. Correlation functions in the different calculational schemes: EIHH(dashed-dotted), FY, 
CRCGV, SVM, HH and NCSM (overlapping curves) 
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